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I. INTRODUCTION 



Quantum mechanics is more rigorous than the classical one, which is the best elabo- 
rated and understood part of physics. Classical mechanics can be thought of as the base 
of quantum mechanics. Lagrangian Mechanics and Hamiltonian mechanics are geometrical 
description of classical mechanics. In physics, quantization is the process of explaining a 
classical understanding of physical phenomena in terms of a newer understanding known as 
quantum mechanics. This is a generalization of the procedure for building quantum mechan- 
ics from classical mechanics. If a physical phenomenon posses Lagrangian or Hamiltonian 
description, the quantization procedure can be easy performed. But for damping classical 
physical phenomena, it is complicated to perform a quantization procedure from the classical 
description to the quantum one. Readers can read the complete history of the important 
ideas in this field in the literature^. Recent works on quantization of damped systems see 
papers^i^-ii^. 

The works of Kochan^^^iii^attracts our attention, because there is close link between 
the works and the classical mechanics. Kochan^ utilized the tools of contact geometry 
to give a new picture of classical mechanics, and then defined a new variational principle 
of least action^^. Based on the classical picture and the variational principle, a quantum 
propagator— for damping systems was derived. 

The ideas of Kochan motivated us to turn our interest to quantum mechanics, because 
we have found a relation between a damping classical mechanical system and conservative 
classical ones, which can be described as a variational principle^. This picture of classical 
mechanics is illustrated by a simple example in Sec. II. Therefore, we attempt to utilize our 
variational principle^ to construct a new quantum propagator. Our approach to define a 
new propagator is explained in Sec. III. In Sec. IV a simple example, where the quantization 
of a particle with friction kx is performed, is reported. In this example, comparison with 
several other approach is made. 
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II. REVIEW ON A NEW VARIATIONAL PRINCIPLE 



A. A Proposition 

To beging, it is necessary to review the new variational principle^^, which is adopted 
to describe a relation between damping classical mechanical system and some conservative 
ones. The relation is represented proposition: 

Proposition II. 1. For any non- conservative classical mechanical system and arbitrary ini- 
tial condition, there exists a conservative system; both systems sharing one and only one 
common phase curve; and the value of the Hamiltonian of the conservative system is equal 
to the sum of the total energy of the non- conservative system on the aforementioned phase 
curve and a constant depending on the initial condition. 

This proposition can be demonstrated by a simple example. Consider a special one- 
dimensional simple mechanical system 

X + = 0, (1) 

where c is a constant. The exact solution of the equation above is 

X = A + T^e"^*, (2) 

where A, rj are constants. Differentiation gives the velocity: 

X = -nrje'^'K (3) 

From the initial condition Xq,Xo, we find A = xq + Xq/k,, rj = —xq/k,. Inverting Eq. ([2]) 
yields 

1 , X — A . . 

t = In (4) 

K r] 

and by substituting into Eq. (j3]), such we have 

X = —n{x — A) (5) 

The dissipative force F in the dissipative system is 

F = Kx. (6) 
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Substituting Eq. ([5]) into Eq. ([6]), the conservative force T is expressed as 

J= = -K^[x - A); 



(7) 



Clearly, the conservative force J-" depends on the initial condition of the dissipative system 
([T]), in other words, an initial condition determines a conservative force. Consequently, a 
new conservative system yields 



x + J^ = Q^x - K^{x - A) = 0. 



(8) 



The stiffness coefficient in this equation must be negative. One can readily verify that the 
particular solution ([2]) of the dissipative system can satisfy the conservative one (|H]). This 
point agrees with Proposition (III. II) . 
The Lagrangian of the system ([8]) is 



1 1 

f— 2,-"- 22 2a 
Lj — — X ~\~ — /t^ X /t^ i\X 

2 2 



(9) 



B. Infinite-dimensional Variational Principle 



The ideas of literatures^""— have motivated us to regard the mechanical system ([T]) as a 
special fluid, which is a collection of fluid particles in phase space. Therefore, first let the 
label of a particle in the phase space be 



a = [Xo,Xo) 



(10) 



the coordinate of a particle in the configuration space 

q = q[a, t) = {x{a, t), x{a, t)} ; 
Po = 1- One can consider L in Eq. ([9]) as a Lagrangian density of the system ([T| 

r f ^-2,^ 2 2 2a 
L = L = —X H K X — K AX 

2 2 

Thus the Lagrangian functional of Eq. ([T]) can be presented as the following: 



L[x, x] = I Cd^a 

' D JD 



-'-•2,-'- 22 2a 

2 2 



d^a, 



where d^ = dxodxo- Thus the action functional can be presented as follows: 



S[q] 



L[x, x]dt 



to 



dt 

to J D 



^■2,^ 2 2 2 a 

-X + -K X — K Aa; 



d^a 



(11) 



(12) 



(13) 



(14) 
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According to Hamiltonian theorem, we have the functional derivative SS/Sq{a,t) = 0, 

6S 6L d 6L 

=^ — - — — = 



6q{a,t) 6q{a,t) dt6q{a,t) 



-x{a,t)-K^{x~A) = 0. (15) 



dq{a,t) dtdq{a,t) 

The equation above imphes that subject to the initial condition a, an associated conservative 
system exists, the control equation of which is Eq. ([8]), the phase curve of which coincides 
with that of the damping system ([T]). In the case of the classical least action principle, 
the stationary curves set is represented as an Euler-Lagrangian differential equation, which 
can be converted into a uniform conservative Newtonian equation. In the case of our least 
action principle, the stationary curves set is represented as an Euler-Lagrangian functional 
equation, which can be converted varied conservative Newtonian equation subject to varied 
initial condition a. 

Furthermore we can define an action density S, such that Eq. ( |T4l) can be rewritten as 



S[q] = I Sd^a = I d^a I Cdt, (16) 

' D Jd J to 



where S = J^^ Cdt. If Eq. ([T]) is conservative, the action density reduces back to the classical 
action, and the Lagrangian density back to a uniform classical Lagrangian, our least action 
principle back to the classical least action principle. 



III. DERIVATION OF A QUANTUM PROPAGATOR 

According to Feynman's theory^, the probability amplitude of the transition of a system 
from the space-time configuration {xa,0) to another space-time configuration {xb,T) is 

/f(xi,r;a:„,0) = i j \V',]expUs\ 

4/polexp{i^Ld«} (17) 

Here the integral is taken over all pathes 7(t) = {x,t) or all pathes in the phase space, 
satisfying the boundary conditions 

x{0)=Xa, x{T) = Xb. (18) 

Kochan employed his variational principle to replace Ldt in Feynman's kernel function with 
a Lepage two-form Q. Motivated by the Kochan's trick^^, encouraged by the sentence from 
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Feynman's thesis^^ too: "the central mathematical concept is the analogue of the action 
in classical mechanics. It is therefore applicable to mechanical systems whose equations of 
motion cannot be put into Hamiltonian form. It is only required that some sort of least 
action principle be available" , we propose a method to generalize the Feynman's probability 
amplitude: 

K(xt,T;x,,0) = ^ J |I>7]exp|is| 

= j^JlV,]e.p{yC.lt], (19) 

where the action S is replaced by the action density S and the Lagrangian L replaced by 
the Lagrangian density C, the classical least action principle is replaced by the variational 
principle in Sec. III. 

In conservative case, the Lagrangian density C reduces back to standard Lagrangian, the 
propagator f|T9l) reduces back to standard Feynman's propagator f|T7|) 



IV. EXAMPLE OF QUANTIZATION 

Consider the quantization of the system (fl]) with the boundary condition (fTSjl . From the 
boundary condition can be derived 

Xh Xa Xb Xa /nr>\ 

We take Feynman's general methocli2. to integrate Eq. (fT9|) . First we construct a short- 
time propagator: 

{Xj^l — Xj) . f 2 I 2\ 2 a /-^i+l -^j ■ 



ir(x,^i,x,;6) = yillexp| 



+ -{x'^^, + x^)-k'A{- 



2e2 4 ^ 



27ri 

where 



^ exp 1^ [ao{x^j+i + x|) - 2boXj+iXj - RqXj+i - SqXj] 



(21) 



Substituting the short-time propagator into 

N-l „+oo ^-1 



/- +00 -IV— -L 

dxk n K{x Xj] e 
., -oo „-_n 



fc=i " j=0 
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we have 



K{xb,T;Xa,0) = lim 

N^oo 



Here 



1 y/' 

2mhe J 



'N-l 



n 



N-l 



X exp ^ XI ["0(^1+1 + ^j) ~ '^boXj+iXj - RoXj+i - SoXj] 



lim Ai 



— lim A2 



— lim Ak 



j=0 
'N-l 



n / 

N-l 

n / 



.i=2 " 
■w-i „ 

n / 



exp -(^1 



exp <; -(^2 



exp <; -0fc 



lim Ajvexp < -0jv 



Ai 



2'Kihe 



N-l 



01 = X [ao(a;|+i + - 26oa;j+iXj] + ai 
Q^i = 00(2^2 + ^1) ~ 26oa;2a;i — R0X2 — SqXi + ao{xl + xl) — 2boXiXQ — RqXi — SqXq 



= 00(^2 + xl) - R0X2 - SqXq + 2ao 
[bo{x2 + xo)^ + {So + Ro)/2f 



Xl 



2an 



Xl - 



2ao 

ho{x2 + xof + {So + Ro)/2 
2an 



(23) 



ho{x2 + Xof+{So + Ro)l2 
2ao 



+ ai{x2 + Xq) — 261X22:0 - R1X2 — SiXo — ^1 



where 

ai = ao - h = T^, (24) 

2ao 2ao 

P ^ , b,{So + Ro) bo{So + Ro) ^ (i^o + Sof 

ni — no-\ bo — oo-\ III — z i^'^J 

zoo 2ao Sao 

Integration is first done with respect to xi, we have 
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N-1 

4>2 = ^2 ho(^j+i + ^j) ~ 2boXj+iXj] + ai{xl + xl) - 261X2X0 - R1X2 - SiXo 

J=2 



N-1 

0fe = ^ [ao(xJ+i + xj) - 2boX-j+iXj] + ak-i{xl + x^) - 2bk-iXkXo - Rk-iXk - Sk-iXg 

j=k 
N-1 

= XI [ao{x^j+i + x'^j) - '^boXj+iXj] + ak, (26) 

j=k+l 

where 

<^k — (^oi^k+i + ^fe) ~ '^boXk+iXk — RoXk^i — SoXk 

+ak-i{xl + Xq) — 2hk^iXkXQ — Rk-iXk — Sk-iXo 
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_ (6oXfc+i + 6fc-iXo) + {Rk-i + Sq)/2 

^k 



[(6oXfc+i + bk-ixp) + (-Rfc-i + '5'o)/2] 
Integration is first done with respect to x^, wc have 



A+i = Ak\l — — — exp i-^^k 



N-1 

= ^ [«o (2^5+1 + x'^j) - 2boXj+iXj] + aoxl^^ + ak-ixl - RoXk+i - Sk-iXo 

j=k+i 

[{boXk+i + bk-ixp) + {Rk-i + SQ)/2f 

N-1 

= [ao{x]^^ + x])-2bQXj+iXj\ 

3=k+l 

+0'kx\j^i + a'fcXo — 26fcXfc+iXo — RkX^+i — SkXo — Q^, (27) 



where 



^0 / ^fc-l 7 bgbk-i , . 

Ofc = oo ^ , dk = \ h = ■ , (28) 

C^O + Qfe-l + 0-k-l Q-0 + 



\2 



P p , (-^fc-1 + '5^o)6o Q Q , (-^fc-1 + So)bk-i (Rk-i + 'S^o) , . 

-n-fc — -n-O H : , -Jfc — -Jfe-l H : , ^^fe ~ 77 T l^^J 

ao + afe_i ao + ajfc_i 4(ao + ajk_i) 

If Cfe = a'^, i.e. 

aU = bl_, + al-bl (30) 
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then the form of (pk+i is same as the form of 0^. Eq. (130|) is tenable as /c = 1, = 2, because 
from Eq. fl25l) can be derived 



2 1,2 



ao 



2a 





2ao 



2 1,2 



Eq. fl27I) and Eq. fl28l) can be proved by mathematical induction for k = 1,2, 3, 



mh /In 

exp ( — — iiAT-i 



ao + aN-2 



h 



1 \ 



2T(ihe 



/ \ I — 1 / I — 1 



!)Ar = aN~l{x\ + a;o) - 26Ar_iXArXo - Rn-iXn - Sn~iXo 



Substituting the equation above into Eq. (12^ . we have 



K{xb,T]Xa,^) = lim 



e-S>0 



1/2 



W-1 

n 

.fc=i 



1 1 

2eao + Ofc-i 



V2 /Ar_i 

exp ^ --r^fe-i 



exp <^ ^ [aN-i{xlj + xl) - 2bN-iXNXo - Rn-iXn - Sn-iXq] \ (31) 



The coefficients a^, bk of Eq. flHT]) can be derived from the following recursion formulas: 



(32) 



(33) 



In order to obtain the result of Eq. fl3T]) with N ^ oo, one must obtain the limit of a^, 
6fc, Sk, fik- Let us consider ak, bk first. As — )■ oo, we have 



2 2 



(34) 



Substituting Eq. (iMj) into Eq. fl32|) . we have 



<^o = ^0 + 2 sinh 2^^J ~ '^^ 
Substituting the equation above into Eq. (133|) . we have 



1 cosh/te 1 / ^0 1 
-r- = —, ^ 1- W 1 + 79 — sinh ne 

bk bk-i bo y bl_^ 



(35) 



Let = 1, Eq. (l35ll becomes 



— = —(cosh Ke + ■>/ 1 — sinh^ Ke) , 
6i 6o bo 



2 cosh Ke sinh 2Ke 



bo sinh 



(36) 
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Let = 2, we have 



1 cosh/te 1 / 6n ,9 sinhS/ce sinhS/te 
— = — ; h i-Wl + T^smh^ = = 

02 bi oq W of Oi sinh 2Ke Oq sinn Ke 



The rest may be deduced by analogy, we have 

1 sinh kne 



Jk-i 



bo sinh Ke 



The equation above can be proved by recursion method. Substituting the equation above 
into Eq. (^^, we have 



Hence 



1 cosh K,e sinh kne 1 
bk-i bo sinh bo 

sinh(A; + l)K,e 
bo sinh ne 



sinh kne 
bo sinh ne 



sinh K,e 



sinh Ke 



2e sinh(/c + l)Ke 



ttk = {bl + al- bl) 



/ fen . , 9 sinh Ke coshffc + l)ft;e 



(37) 
(38) 



1 / coshfc/te 
ao + cifc-i = — coshne + smh /te — — — 
2e \ smh fc/te 

1 sinh(A; + l)K,e 



2e 



sm Ke 



(39) 



Substituting the equation above into Eq. (1291) and let = 1, we have 

ao + ao ^ V "^osh / 



1 



ao + ao 2 \^ ' cosh ne 

Let k = 2 and substitute Eq. (l39ll into Eq. (l29l) . we have 



i?2 

^2 



+ S'o)6o /«'Ae 



ao + ai 



-Ri + 5o , 

Si H ^ bl 

ao + ai 



1 + 



1 + 



1 (2coshKe + 1) sinh2Ke 



cosh Ke 
1 

cosh Ke 



1 + 



sinh Sue 
(2 cosh ne + 1) sinh Ke 
sinh 3Ke 
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For k = 3, we have 



ao + 0-2 ^ 

C _ C I -^2 + 'S'o 

•Js — '-'2 H ; O2 

K^Ae r 1 



1 + 



1 2 cosh ne sinh Sue + (2 cosh /te + 1) sinh 2/te 



cosh Ke 



sinh 4Ke 



cosh /te 



1 + (2 cosh(Ke + 1) 



sinh K,e 
sinh Sue 



2 cosh K,e sinh Skc + (2 cosh ne + 1) sinh 2K,e sinh /te 



sinh 3fi;e 



sinh 4fi;e 



We can suppose 



Rk 



1 + 



1 X]j=i 2 cosh /te sinh j/te 
cosh/te sinh(/c + l)Ke 



(40) 



1 



sinh Ke 



cosh Ke sinh(j + l)/te 



'Yl,l=i 2 cosh Ke sinh lut 
sinh j/te 



(41) 



Substituting Eq. fHOj) into the last term in Eq. fl29|l . we have 



^3^2g3 gjj^]^ 

4 smh.{k + l)Ke 



1 X]y=i 2 cosh Ke sinh jKe 



2 + 



cosh 



sinh /cKe 



(42) 



Therefore, we have the hmit of Rk, Sk, 



lim R - ^'Ae^"^-2e"^ + l _ ^^^^^ e"^ - 1 
Af-i>oo ^ 2 Ke'^^ sinh kT e'^'^ + 1 



(43) 



hm Sn = kA 



(44) 



N-l 



hm Vfi 

Af-i>oo 



k-1 



k=l 



(45) 



Substituting Eq. §^ and Eq. ([37]), Eq. ([391), Eq. (03]), Eq. (|l3]), Eq. (|45]) into Eq. ([23]), 
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we have 



K{xb,T;Xa,0) = lim 



X exp 



1/2 



'N-l 

n 

.fc=i 



1 



2eao + ak-i 



1/2 



exp 



''N-l . 

5Z 



h 



lim 

N^oo 



2iTihe 



1/2 



n 



sinh Ke 



exp^- 



ysinh(A; + l)Ke 
sinh Ke cosh A^Ke 



1/2 



N-1 



exp hm > i7 



fc-1 



fc=l 



2 2\ l Si'^h '^^ 
(x^ + Xq) . , — ; — 7^; X]\fXo 



2e sinhA^Ke^ ' e sinh(/c + l)Ke 

— hm RnXn — hm S'atXq ? 



27ri/isinh kT 



1/2 



exp < — 



h i2 sinh kT 



-kA- 



[Xb + Xa 



Substituting Eq. (fTSjl into the Equation above, we have 

1/2 



K{xb,T;Xa,0) 



K 



27cih sinh kT 



exp < — 



h i2 sinh kT 



[cosh kT{xI + xl) — 2xbXa\ (46) 

(47) 

[cosh kT{xI + xl) — 2XbXa\ 



-kT 



K I Xa 

2mh sinh kT 



<^ b a 1 ^ 1 



- 1 / e'^^ + 1 



1/2 



(a;fe + Xa) 



exp^- 



2 tanh kT 



{x1 + x^) 



KXbXa 

sinh kT 



-K X, 



- 1 



tanh— (xf, + Xa) 



(48) 



It is evident that the propagator above tends to the free Schrodinger propagator in the hmit 
K ^ 0. 

Finahy we perform an analysis of the time evolution in terms of the quantum propagator 
( 148|) -Suppose that at the initial time we have a Gaussian wave-packet 

ipQ^q) oc exp ( -6*0?^ + ^v^q 



h 

, which describes a unit mass particle localized in a neighborhood of the point g = with 
an initial velocity vq. At a later time t, the system under consideration will be characterized 
by the propagated wave-packet distribution, i.e. 



i)T{x) 



K{x, T; q, 0)tpo{q)dq oc exp { -6*1 (x - (x))^ + -6'2g^ 



(49) 
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The direct calculation shows that the Wave-function i/^rix) still remains Gaussian. The 
mean value of position is 

(.) = .0^^^ (50) 

K 

The mean value of velocity of the wave packet 



, , 2k tanh (^) , , 

{v)= ^_,^^lU ^)+Vo (51) 



The coefficient Oi in Eq. fl49|) 

tanh^(K*T) 



9i 



"0 



n tanh ^ — n 



tanh ^ 



(52) 



Let us compare our result with the results obtained by Kochan's approach^>^ and 
Caldirola-Kanai's methocl^, Das's methocl^. The results obtained by Kochan's (abbr. 
Koch) approach^ are 

tanh 



{x)Koch = '2vo- 



2 



K 



{v)Koch = - 2^(^) 



^ _ 16^12 tanh^{0.5KT) 



_ ■ k(3-c-''^) 
^ 4fi(l-e'=T) 



The results obtained by Caldirola-Kanai (abbr. CK) method are 

K 

{v)cK = Vo 

4ri2(l-e-«J')2 



"0 ^2^(1-6-"'^') 

The results obtained by Das's (abbr. DGST) method are 

(a:WT = ^(l-e-«^) 

\v)dgst — ^oe 



The comparison among these results above is illustrated in Fig. [H where q = Om, vq = 
5ms~^, K, = 0.6s^^ . The curve obtained by our method is labeled with 'LG'. It is clear that 
the expectation value of velocity obtained by our method and Kochan's method look very 
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FIG. 1. Comparison of the mean values of position and velocity obtained by various methods 

similar. The expectation value of position obtained by our method is quite different from 
that obtained by Kochan's method, but with time evolution, (x) obtained by our method 
approaches {x)ck- The reason is 



T—>-oo T— >oo K T— >oo K 
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Because our averaged velocity (v) and Kochan's averaged velocity posses a similar advantage 
that they decrease with respect to time translations, also possess similar restriction that our 
averaged velocity is reliable in the time interval [0,ln(-\/2 + 1)//^] and Kochan's reliable 
interval is [0,ln3/fi;]. Exactly at the moment ln(v^ + 1)/'^, our averaged velocity becomes 
zero and going over, it would produce a negative (nonphysical) mean velocity. 



V. CONCLUSION 

Kochanii have compared his results with other results obtained by the method of 
Caldirola-Kanai and other methods. Because the aforementioned similarity, our approach 
possesses the advantage of Kochan's approach, but the position obtained by our approach 
likes that obtained by CK-method. Since both Kochan's approach and our approach are 
based ultimately on the classical history, this similarity exists between the approaches. For 
quantization of more complicated damping mechanical systems, one can utilize numeri- 
cal method to obtain classical trajectory, then construct approximation of the Lagrangian 
density. 
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